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We review briefly the current status of tliermal pliase transitions witliin tlie 
Standard Model and its simplest extensions. We start with an update on QCD 
thermodynamics, then discuss the electroweak phase transition, particularly 
in supersymmetric extensions of the Standard Model, and end with a few 
remarks on the cosmological constraints that thermal phase transitions might 
impose on even higher scale particle physics. 
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1 Introduction 



As is well known, the microscopic energy scales of quantum mechanics and the macro- 
scopic properties of our present Universe are intimately connected. For instance, the 
0{eV) energy scale of atomic physics manifests itself through the existence of the cosmic 
microwave background radiation, and the 0{MeV) scale of nuclear physics through the 
primordial origin of light element abundances. The connection might of course extend 
even much further on: the small fluctuations observed in the microwave background could 
have been produced during an early period of inflation, which could well be a manifesta- 
tion of particle physics at, say, the 0(10^^ GeV) scale of grand unification. 

With this background, it is natural to expect that the main scales of experimentally 
accessible particle physics should also leave their marks in cosmology. Could not the 
C(GeV) scale of QCD, and the C(TeV) scale of the electroweak (EW) theory, be related 
to such cosmological remnants as extragalactic magnetic fields, or baryon asymmetry? 

Taking a closer look, it turns out that quite non-trivial conditions have to be met for 
the latter connections to exist. Among the biggest challenges are: 

1. Generic QCD and EW interactions are so strong that it is difficult to deviate 
from thermal equilibrium, which is a necessary requirement for any cosmological remnant 
to emerge. Denoting by n a particle density, a a cross section, v an average velocity, 
a = 5'^/(47r) the coupling, T the temperature, and mpi ~ 10^^ GeV the Planck mass, 
interactions fall out of equilibrium provided that 

1 1 mpi 

T = :^t = -^, (1) 

nav a'^i I ^ 

which leads to T > a'^m-pi ~ 10^^ GeV. Thus, they can fall out of equilibrium only at 
temperatures where anyway the QCD and EW theories merge to a grand unified theory. 

There is one major exception to the conclusion just drawn: a theory possessing a first 
order phase transition falls out of equilibrium, even if microscopic interactions are strong 
enough to be in local thermal equilibrium above and below the transition point. 

2. Spatial fiuctuations (relevant to remnants such as magnetic fields) are produced 
only on very small length scales. Indeed, the horizon of the moment when T ~ 1 GeV 
corresponds today to about 1 light year, and that of T ~ 1 TeV to about 1 astronomical 
unit. Fluctuations can effectively only be produced on scales smaller than these, which 
leads to miniscule numbers with respect to intergalactic distances. 

Again, there is one conceivable way of avoiding the problem: magnetohydrodynamic 
evolution is very non-linear and could potentially transfer energy to large length scales 
more rapidly than comoving expansion 

The subject of this talk is the analysis of the first of the problems mentioned, the 
existence of first order phase transitions in QCD and in various versions of the EW 
theory. The actual generation of cosmological remnants has been treated in other talks 
at this conference H, H], and will only be touched upon very briefiy here. 
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2 QCD thermodynamics 



The theory of strong interactions, QCD, is expected to undergo a phase transition at a 
temperature of the order of A-^. The transition is said to be related to deconfinement and 
chiral symmetry restoration. However, neither of these are rigorous concepts for physical 
quark masses, and therefore there might just as well be a smooth gradual change of the 
properties of the system, instead of an actual singularity. 

In any case, a smooth change is not what one naively expects. Indeed, counting 
the "free" degrees of freedom in the "confined" and "deconfined" phases, one finds a 
considerable change, suggesting perhaps a strong transition: 



However, in practice, particles are not free but strongly interacting. This leads to the 
fact that the properties of the QCD phase transition can, on a quantitative level, only be 
studied systematically with four-dimensional (4d) finite temperature lattice simulations. 

It turns out, furthermore, that such lattice simulations are very demanding. The rea- 
son is that the precise characteristics of the transition, such as its order, depend strongly 
on the symmetries of the system, which are in turn determined by the quark masses (for 
a review, see Q). But chiral quarks, as light as they are in Nature, are difficult to fit on 
the lattices available in practice. Thus it is believed that the properties of the system do 
change at a temperature Tc ~ 170 MeV [^], but whether the change is smooth or, in the 
large volume limit, discontinuous, remains open. 

What is known much better is the behaviour of various thermodynamical quantities, 
such as the pressure, at temperatures above the critical, T > Tc. Indeed, there the inclu- 
sion of quarks does not appear to change the result qualitatively from the idealised 
case of pure SU(3) 0. The general pattern, illustrated in Fig. |], is that the pressure is 
small at T ~ Tc, rises rapidly at T ~ (1...2)Tc, and then levels off, approaching the ideal 
gas limit very slowly, with a deviation of 10. ..15% for a long while. 

This noticeable deviation from ideal gas thermodynamics has some implications. Ac- 
cording to the Einstein equations, the temperature in the Universe evolves as 



Here s = p'{T),e = T^{p(T)/Ty . Using an equation of state as in Fig. |I|, the Universe 
cools slower than a free gas 0, and the sound velocity = p'(T)/e'(T), characterising 
hydrodynamic fiuctuations, dips around the transition point. These facts do result in a 
non-trivial fiuctuation spectrum, even if long-lasting consequences might not remain. 
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Figure 1: The pressure of pure SU(3) gauge theory, compared with the free Stefan- 
Boltzmann value po = (7r^/45)(A^^ — 1)T^, as a function of T/A^g, where Ams is the scale 
parameter. The transition takes place at T ^ A^g- The 4d lattice results are from , and 
the curve labelled perturbation theory + 3d lattice, from 0. The results are compared 
with the "bag" equation of state, pba_g{T) = Po(T) — pq(Tc). 



3 EW phase transition 



At temperatures of the order of 

weak coupling constant, the electroweak gauge symmetry gets restored 



mff/g, where rriH is the Higgs mass and g is the 

TD|. This may 



have led to an important consequence, the existence of a matter-antimatter (baryon) 
asymmetry |TT|] (for reviews, see |T^, |T^, 0). Again, a deviation from equilibrium through 
a first order phase transition is a necessary requirement. More quantitatively, a 1-loop 
saddle point computation [14| as well as a non-perturbative evaluation |15| show that the 
discontinuity in the Higgs expectation value across the phase transition, Av/T (in, say, 
the Landau gauge) should exceed unity, Av/T> 1.0. 

In contrast to the QCD case, many properties of the EW transition can be addressed 
in perturbation theory. This is simply because the Higgs mechanism itself is perturbative. 
Eventually perturbation theory breaks down, though: at finite temperatures the largest 



loop expansion parameter can be estimated to be |16 



nm 



(5) 



where m is some mass scale. Thus light excitations, m < g'^T, always present at the phase 
transition point, lead to an infrared problem . 

The expansions parameters related to heavy excitations, on the other hand, such as 
the non-zero Matsubara modes of finite temperature field theory with m ~ vrT, are small. 
This allows for an essential simplification of the non-perturbative treatment: one can 
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Figure 2: The phase diagram of the physical Standard Model. The blobs indicate simu- 
lation points (with error bars) ||2^, |25|, and the solid curve is a fit through them. 



integrate out massive modes perturbatively and study only the dynamics of the light 
modes non-perturbatively [|18|. In the first step, the original 4d theory reduces to a three- 



dimensional (3d) effective one, in a process called dimensional reduction [|T^ . For studying 



the EW phase transition in a weakly coupled theory {rriH < 250 GeV), this approach works 
with a practical accuracy at the per cent level, from the point of view of both dimensional 

PI, 



reduction ||2^-||2^ and numerical simulations 

In the case of the Standard Model and many of its extensions |^, the only non- 
perturbative infrared modes are a Higgs doublet and the spatial SU(2) and U(l) gauge 
fields, with field strength tensors Fij,Bij and gauge couplings g3,g'^- The Lagrangian is 



3d 



^TrF^. 
2 



4 + 



+ A3 



(6) 



All knowledge about the physical zero temperature parameters, as well as about the 
temperature, is encoded in the expressions for the effective couplings ml, XziQhd'i- 

The properties of the phase transition now depend on the effective couplings in Eq. (P). 
In particular, since the gauge couplings are fixed and the mass parameter m\ is to be tuned 
to the phase transition point, such properties are determined by the single dimensionless 
ratio X^l g1 To get a feeling of this dependence, let us first apply 1-loop perturbation 
ori 

discontinuity 



theory. Ignoring the tiny corrections from g'^ j g\, we find a first order transition with the 



gT 



9l 



1 

8^ A3 



^gi^ 



(7) 



Thus, the transition weakens for large Xs/g^- In reality, the line of first order transitions 
even ends completely, at Xs/g^ = 0.0983(15) 
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3.1 Standard Model 



To be now more specific, let us consider the Standard Model. Then 



+ O 

gl 8m^ V(47r)2m' 



Applying this (with the actual 1-loop corrections properly included) to the lattice re- 
sults of @ |2| leads to Fig. |. We find that the endpoint location \^/gl = 0.0983(15) 
corresponds in physical units to mH,c = 72.3(7) GeV, = 109.2(8) GeV [^ . 

As pointed out above, for baryogenesis we would need a first order phase transition, 
and even a strong one, Av/T > 1. We have learned that a first order transition only exists 
for rriH < 72 GeV. Furthermore, it has in fact Av/T < 1.0 down to mn ~ 10 GeV El 



Thus, we observe that experimentally allowed Higgs masses niH ^115 GeV ||2^ are very 
far from allowing for electroweak baryogenesis. 

It appears that primordial magnetic fields present at the time of the EW transition can 
strengthen the transition quite significantly, but not enough to change the conclusions i28[| . 



though they are associated with other intriguing phenomena |29 . 



3.2 Supersymmetric extensions of the Standard Model 

How should the Higgs sector be modified in order to change the pattern above? As 
the strength of the transition is determined by the scalar self- coupling, cf. Eq. (|^), we 
apparently need some new degree of freedom, which can decrease the effective A3 by 
0(100%). It turns out that to get such a large correction, we need a bosonic Matsubara 
zero mode with an expansion parameter as in Eq. (|^). A simple perturbative 1-loop 
computation shows that at least the sign of such loop corrections is the correct one, 
negative. But to have an effect of (9(100%), we need e ~ 1, so that m ~ g'^T/n. That is, 
the degree of freedom should itself be non-perturbative! 

In the MSSM, natural candidates for such degrees of freedom are squarks and sleptons. 
Then there are thermal corrections in their effective mass parameters, appearing as mg ~ 
"^4d + 6'^^^- Thus, to get a total outcome of order (gf^T)^, a cancellation must take place, 
which requires a negative mass parameter m1^. At zero temperature, the physical mass 
is then roughly ^«phys~"^4d + "^top~'^top- order for such a relatively light degree 
of freedom not to have shown up so far in precision observables, it should be an SU(2) 
singlet. Since it should also couple strongly to the Higgs, we should choose stops. 

Now, the stops come left- and right-handed, tL,^^ (these states can also mix, but 
for lack of space we ignore this here). The requirement of having a small violation of 
the electroweak precision observables means that the weakly interacting left-handed one 
cannot be light, m^^ 3> mtop- This also increases the Higgs mass (see, e.g., pO|), 

ml ~ m| cos^2/3 + |44^ 1^ (9) 
87r^ mf^ mtop 
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Figure 3: Examples of values (Higgs mass mn, mass parameter rriA), leading to a large 
Av/T according to 2-loop perturbation theory, for an optimal choice of m^^ (from ||35|| ). 
The region indicated with arrows (widening for increasing m^^ |33, 35 1) appears to be 
within experimental constraints, even for the most stringent "no mixing" scenario 



On the contrary, the SU(2) singlet stop iji can be "light", and serve as the desired new 
degree of freedom [^]. (Then, of course, we lose half of the correction to rriH in Eq. (^, 
but this is the price to pay.) It should not be so light that the stop direction gets broken 
before the phase transition, though, because then one cannot get back to the EW minimum 
afterwards |^ . Another example of a viable light scalar degree of freedom is the complex 
gauge singlet of the NMSSM 13 . 

Under the conditions described above, resummed 2-loop perturbation theory indicates 



that the EW phase transition can indeed be strong enough for baryogenesis |34]. To check 
the reliability of perturbation theory, a 3d effective theory (a generalisation of Eq. (^) has 
again been constructed and studied with simulations with the encouraging outcome 
that, in fact, for a light right-handed stop and heavy Higgs, the 2-loop estimates are 
reliable and even somewhat conservative. This is in strong contrast to the case of the 
Standard Model at realistic Higgs masses. 

Next, we must ask whether this parameter region is indeed in agreement with all ex- 
perimental data. The most important constraint comes from the lower bound on the Higgs 
mass There is a parameter in the MSSM, rriA, which determines whether the Higgs 
sector resembles that in the Standard Model (m^ > 120 GeV) or not (m^ < 120 GeV). 
In the latter case, the experimental lower bound is relaxed |2^, and since the transition 
needs not always get significantly weaker (see Fig. this case is acceptable. If is 
larger, then the left-handed stop should be quite heavy, m^^ > 2 TeV, in order to increase 
the Higgs mass towards the Standard Model value > 115 GeV (cf. Eq. (|)). Allowing for 
significant mixing in the stop mass matrix relaxes the Higgs mass bounds quite consider- 
ably although it also tends to weaken the transition somewhat [0, p5|] . 
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4 Phase transitions at higher scales? 



Let us end by considering hypothetical phase transitions at scales even higher than the 
EW one, possibly related to unification. To discuss this self-consistently, it must be 
assumed that infiation and unification are not related, but infiation takes place earlier 
on, maybe at the Planck scale. Then, in clear contrast to the cases considered so far, 
thermal phase transitions rather generically do tend to produce cosmological remnants, 
whose non-observation allows to place constraints on possible unification models, or on 
cosmology. We are here referring to topological defects, such as domain walls, cosmic 
strings, and monopoles [36|. As an example, let us make some more specific comments 
on the first ones. 

The reason why there is a strong constraint is that the energy density 5p related to 
domain walls decays much more slowly than radiation. It would lead to 5p/ p > 10~^ 
at photon decoupling, if the domain wall surface energy density is cr > (1 MeV)^. This 
exceeds the fiuctuations observed in the cosmic microwave background. Thus any new 
theories leading to thermal domain wall production are excluded [P^! 

As timely examples, let us consider models with compact extra dimensions. In case 
there are gauge fields living in the bulk ||3^ but no fundamental matter, one generically 
gets Z{Nc) domain walls [^]. Thus, a class of such models may be constrained by the 
thermal phase transitions they would undergo in cosmology. It should also be mentioned 
that the winding of a brane around the extra dimension may lead to various types of 
cosmic strings and monopoles EO], potentially resulting in other constraints. 



5 Conclusions 

To summarise, the QCD phase transition is probably very weak, if there at all. However, 
the system is strongly interacting, and its thermodynamics deviates significantly from 
that of an ideal gas, which leads to a lengthy period of non-standard expansion. 

In EW theories, there is in general no phase transition at all for realistic Higgs masses, 
unless there is also another light scalar degree of freedom, which plays an essential role 
in phase transition dynamics. For instance, the transition can be strong enough for 
baryogenesis in the MSSM if there is one very light {tTL^^ < rritop) and one rather heavy 
stop. Either the heavy stop should be really heavy, m^^ ^lOmtop, or the Higgs sector 
should contain at least two independent light particles, niA < 120 GeV, in order not to 
violate experimental constraints. There is much more freedom in the NMSSM. 

While little concrete can be said about theories of unification scale physics, one can 
at least note that the possible overproduction of topological defects in phase transitions 
can place some constraints on model building and cosmology (on top of many other 
constraints, of course). A classic example is the monopole problem, but one can also 
arrive at a domain wall problem in some models with extra dimensions. Such problems 
can be avoided if the reheating temperature after infiation is below the unification scale. 
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Finally, let us stress that the topic of this talk has been phase transitions taking place 
in a system in local thermodynamical equilibrium. In case of low-scale inflation ending 
at a reheating temperature T ~ 0{TeV), it is also natural to consider non-thermal phase 
transitions (see, e.g., [0), which lead to many new physics possibilities. 
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